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In- this  vnpq.rt,  wg  ^review^a  method  for  computing  the  aerodynamic 
influence  coefficients  (AICs)  for  surfaces  with  supersonic  leading  edges.  The 
method  is  based  on  the  two-dimensional  second-order  potential  solution  of 
Van  Dyke.  The  strip  oscillatory  coefficients  are  obtained  from  the  extension 
of  Van  Dyke's  work  by  Rodden  and  Revell  to  include  the  effects  of  sweep  and 
finite  span  in  addition  to  the  effects  of  thickness,  leading  to  a  method  that  will 
avoid  the  unconservatism  of  linearized  theory  and  wilL  be  applicable  at  Mach 
numbers  below  the  lower  Mach  number  limit  of  piston  theory. 


The  influence  coefficients  relate  the  aerodynamic  forces  to  the  surface 
deflections^through  .the  following  definitions.  In  the  oscillatory  case, 

,/  |f|  =  pwbr  s  [Ch]  |hf 

/ 

/ 

and  i yi  the  steady  case, 

/ 

/  |Fs|  =  (l/Z)pV2(S/c)[Chs]  \h\ 


•The  Aerospace  IBM  7090  Computer  Program  No.  HM10  provides  the  AICs  in 

'printed  and  optional  punched-card  output  formats.  The  theoretical  formula- 

i 

I  tion  is  limited  to  Mach  numbers  normal  to  the  leading  edge  greater  than  1.25. 
^The  program  capacity  is  25  surface  strips,  25  values  of  Mach  number,  and  50 
values  of  reduced  velocity. 
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SYMBOLS 


Coefficients  in  series  expansion  for  oscillatory 
aerodynamic  coefficients 

Coefficients  in  series  expansion  of  pressure 
coefficients 

Local  semichord 

Reference  semichord 

Element  of  oscillatory  aerodynamic  influence 
coefficient  matrix 

Element  of  steady  aerodynamic  influence  coeffi¬ 
cient  matrix 

lyocal  chord 

! 

Mean  aerodynamic  chord 
Floot  chord 
Tip  chord 

Distance  between  forward  and  aft  control  points 

Control  point  force 

Finite  span  correction  factors 

Airfoil  semithickness 

Vertical  deflection 

Thickness  integrals 

Local  reduced  frequency,  k  =  ub/V 

Reference  reduced  frequency 

Oscillatory  leading  edge  lift  coefficients 

Lift  referred  to  leading  edge  motion 

Freestream  Mach  number 


v 


M 

Normal  Mach  number,  M  =  M  cos  A 

Mu  .  M 

Oscillatory  leading  edge  moment  coefficients 

o  o 

M 

o 

Moment  about  leading  edge 

N 

=  [(V+l)/2]  (M/p)2 

N  • 

=  [(V+l)/2]  (M/p)2 

r 

=  Tte /t 

S 

Wing  area 

s 

Wing  semispan 

V 

Free  stream  velocity 

V/b  w 
r 

Reference  reduced  velocity,  V/br<*>  =  l/hr 

x,  £ 

Streamwise  coordinates 

X 

Chordwise  location  of  point  of  maximum  thickness 

m 

a 

Angle  of  attack 

P 

=  (M2  -  1)1/2 

P 

=  (M2  .  1)1/2 

y 

Specific  heat  ratio  of  air  (y  =  1.400) 

Ay 

Strip  width 

A 

Leading  edge  sweep  angle 

X. 

Trailing  edge  sweep  angle 

p 

Free  stream  density 

T 

Airfoil  maximum  thickness  ratio 

Tte 

4 

Airfoil  trailing  edge  thickness  ratio 

<*> 

Circular  frequency 

[  J 

Square  matrix 

1  1 

Column  matrix 

vi 

SECTION  I 


FORMULATION  OF  PROBLEM 


A.  Introduction 

For  aeroelastic  analyses  of  surfaces  with  supersonic  leading  edges,  a 
strip  theory  is  desirable  that  will  avoid  the  uncone ervatism  of  linearized 
theory  and  will  be  applicable  at  lower  Mach  numbers  than  piston  theory.  The 
second-order  theory  of  Van  Dyke*  offers  such  a  possibility,  and  it  can  be 
shown^  that  both  the  linearized  theory  and  the  piston  theory  are  special  cases 
of  this  solution. 

Van  Dyke's  analysis  is  the  result  of  both  an  iteration  and  a  frequency 
expansion  in  the  solution  of  the  nonlinear,  unsteady  potential  equation.  The 
iteration  has  been  carried  through  to  obtain  the  second-order  velocity  poten¬ 
tial  (first-order  in  thickness),  and  the  frequency  expansion  has  been  extended 
to  the  cubic  term.  The  result  has  been  presented  in  a  series  form  for  the 
local  pressure  coefficient  on  the  airfoil;  no  control  surface  has  been  con- 
side red. 

The  present  formulation  makes  use  of  the  pressure  coefficient  as  in¬ 
tegrated  by  Rodden  and  Revell,  c  corrected  for  sweep  and  finite  span.  The 
AlCs  for  strip  theory  have  been  discussed  in  Ref.  3  considering  quarter-chord 
coefficients;  a  more  convenient  derivation  utilizing  leading  edge  coefficients 
is  given  in  the  present  study.  The  computational  aspects  of  this  report  are 
an  extension  of  work  previously  discussed  in  Ref.  4. 

B.  Sign  Convention 

The  flutter  sign  convention  is  used  in  the  oscillatory  case:  forces  and 
deflections  are  positive  down;  rotations  are  positive  with  leading  edge  up. 

The  aerodynamic  sign  convention  is  used  in  the  steady  case:  forces  and  de¬ 
flections  are  positive  up,  rotations  are  positive  with  leading  edge  up. 
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C .  Derivation  of  Equations 

The  derivation  is  based  on  the  defining  equation  for  the  matrix  of 
oscillatory  AICs 


jFf  =  pW2  b/s  [Ch]  jh} 


The  steady  AICs  are  defined  by 


|Fs}  =  (1/2)  pV2(S/c)  [Ch8]  {h} 


and  will  be  obtained  from  the  limit  of  the  oscillatory  case 

[ChB]  =  kU“0  2kr2  (8a/S)  ^h1 


We  shall  derive  the  aerodynamic  matrix  for  only  a  single  strip  since  strip 
theory  leads  to  a  partitioned  form  for  the  entire  surface;  e.  g.  , 


in  the  case  of  a  three-strip  wing  where  the  first  null  partition  is  reserved 
for  control  points  whose  aerodynamic  forces  can  be  neglected  (e.g.  ,  external 
stores).  Each  surface  partition  is  of  the  order  2x2  since  we  consider  only 


2 


the  degrees  of  freedom  in  pitching  and  plunging:  the  Van  Dyke  solution  has 
not  yet  been  extended  to  the  control  surface  degree  of  freedom.  The  given 
and  equivalent  force  systems  and  geometry  are  shown  in  Fig.  1.  The  forward 
control  point  in  the  equivalent  system  has  been  arbitrarily  placed  in  the 
quartdr  chord  location. 

The  force  equivalence  ia  given  by 


r 1  1  l 

!F>! 

(  L 

1  o 

Lb/2  (d+b/2)  J 

~  ]  M 
'  o ; 

from  which 

fFl)  R  1+b/ 2d)  -1/dl  /LI 

(F2)  L  -b/2d  1/d  J  (Mq) 

The  leading  edge  oscillatory  coefficients  are  defined  by 


fci  ■  4p“2b2Ay  [o  3 

[Lh0 

K0 

L.ol  jho 
Ma<)J  (ba 

The  geometrical  equivalence  is 

JhJ  ["(1+b/  2d) 

-b/  2d  ”| 

hi 

|ba)  =  |_  -b/d 

b/d  J 

ki 

3 


Fig.  1.  Original  (a)  and  Equivalent  (b) 
Force  Systems  and  Geometry. 
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Cascading  the  above  solutions  we  find 


(F  )  ["(  l+b/ 2d)  -b/d] 

(f‘|  =  4»“  b  ^  [  ,b/2d  b/dJ 


K 

L 

CL 

"( l+b/  2d) 

-b/  2d" 

o 

O 

Mh 

M 

a 

-b/d 

b/d 

o  oJ 


Comparing  this  result  to  the  definition  of  the  influence  coefficient  matrix, 
we  find 


4(b/br)2 

(Ay/s) 

f(l+b/2d)  -b/d] 

L  -b/2d  b/dj 

[Lh 

o 

L 

Qo 

“(l+b/ 2d)  -b/  2d" 

X 

“h 

L  o 

M 

Q 

oJ 

-b/d  b/d 

We  note  the  pre-  and  post-multipliers  to  be  mutual  transposes.  Once  the 
partitions  have  been  obtained  for  each  strip,  the  total  matrix  is  formed  as 
indicated  above'.' 

The  steady  influence  coefficients  follow  as  a  limiting  case  of  the  oscil¬ 
latory  matrix.  From  the  definitions  of  both  oscillatory  and  steady  matrices 
we  find  (since  in  both  sign  conventions  the  forces  and  deflections  are  in  the 
same  direction) 


[ctJ  -  k“s!o 2k/  <«/s>  [ChJ  • 

r 
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We  next  consider  the  oscillatory  coefficients  necessary  in  the  influence 
coefficient  formulation.  Space  permits  only  a  summary  of  the  results  given 
by  Ref.  2.  The  strip  leading  edge  oscillatory  coefficients  are  given  by 


Lh  =  A2f2  +  +  A4f3k) 

o 

L  =  (B.f./k2  +  B-f-J  +  i(B,f,/k  +  B .f.k) 
a  11  Z  j  3  L  44 

o 

Mh  =  C2f3  +  i(G3£2/k  +  C4f4k) 
o 

Mq  =  (Djf2  /k2  +  D2f4)  +  i(D3f3/k  4  D4f5k) 


where  A  ,  B  ,  C  ,  and  D  depend  on  the  Mach  number  and  airfoil  thick- 
n  n  n  n  r 

ness  integrals,  the  fR  are  finite  span  correction  factors,  and  the  reduced 
frequency  k  is  based  on  the  local  chord. 

The  basic  expressions  involved  in  A  ,  B  ,  C  ,  and  D  can  be  outlined, 
r  n  n  n  n 

but  the  reader  is  referred  to  Ref.  2  for  the  finite  span  correction  factors  f^ 

(it  may  be  noted  that  the  f  =  1  for  all  strips  inboard  of  the  tip  Mach  cone, 
where  as  the  f  <  1  in  the  region  of  the  tip  Mach  cone).  The  basic  expres¬ 
sions  are 


II 

CM 

< 

-(1/8)  (a22/2 

+  Jla24  H 

I2a26) 

A3  = 

-(1/8)  (a12  + 

Ila14> 

> 

4*. 

II 

-(1/8)  (a32/3 

+  Kia34 

+  J2a36  +  I3a38) 

B1  = 

(1/4)  (aQl  +  I 

La02) 

Bz  = 

(1/4)  (a21/3  + 

Kla23  + 

JZaZ5  +  I3aZ7^ 

B3  = 

(1/4)  (an/2  + 

Jlal 3  + 

X2a15* 

B4  = 

(1/4)  (a31/4  + 

Lla33  + 

KZa35  +  J3a37  +  X4a39^ 
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where 


C2  =  -0/4)  (a22/3  +  J2a24  +  I3a26) 

C3  =  -(1/4)  (a12/2  +  I2a14) 

C4  =  -(i/4)  (a32/4  +  K2a34  +  J3a36  +  !4a38^ 
D1  =  d/2)  (aoi/2  +  I2a02> 

°2  =  ^a21^4  +  K2a23  +  J3a25  +  I4a27^ 

D3  =  (1/2)  (an/3  +  J2an  +  I3a15) 

D4  =  2 ^  (a31  ^  +  lj2a33  +  K 3a 3 5  +  J4a37  + 

a01  =  "  cos  A 
a02  =  “  (4/P2)  ’  2) 

au  =  (8/^3)  (2  -  M2)  cos  A 

a12  =  (®/^  cos  A 
an  =  (16MZ/p4)  (N  1) 

au  =  (8/p2)  (M2N  -  2) 
a15  =  (8/p4)  (2  -  M2)  (M2N  -  1) 
a21  =  (4/p5)  (M2  +  2)  cosA 
a22  =  O^/P3)  cos  A 

a23  =  (8M2/P6)  [  3  (3M2  -  2)  N  -  2  (5M2  -  3)] 


a24  =  (8M2/p4)  (4N  -  5) 

a25  =  (4M2/p6)  [  <16  -  7M2)  N  +  4  (2M2  -  3)] 

a26  =  (16M2/p4)  <N  -  1) 

a2?  =  (4M2/(56)  [<M2  +  2)  N  -  4] 

a31  =  "  (8M2/3(3^)  (M2  +  4)  cosA 

a32  =  -  (24M2/p5)  cos  A 

a33  =  (16M2/p8)  [  (1 7M4  -  10M2  -  4)  N  -  (5M2  -  2)  (4M2  -  1)] 
=  (16MZ/06)  [  6M2  -  (5M2  -  2)  N] 
a35  =  <16M2/p8)  [7M2  (2M2  -  1)  -  (12M4  -  3M2  -  4)  N] 

a36  =  (8M2/P6)  [2(M2  +  1)  -  (M2  +  8)  N] 

a3?  =  (8M2/p8)  [<3M2  +  2)  -  (3M2  +  4)  N] 

a3g  =  (8M2/(36)  [2(M2  +  1)  -  3M2N] 

a39  =  (8M2/3p8)  [<5M2  +  2)  -  M2N  (M2  +  4)] 
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and  the  thickness  integrals  are  defined  by 


-2b 


I  =  (l/2b)n  f  xn"1g'dx,  n  =  1(1)5 

FI  ^  ft 


0 


J  = 


n 


•n+1  fb  /-1 


(l/2b)1Ti  J  x“  *g  dx,  n  =  1(1)4 
0 


-2b 


K 


n 


=  (l/2b)n+2/  xn  *dxj  gd£ ,  n  =  1(1)3 


.2b 


L  = 
n 


(1  / 2b )n+3  /  xn  ^ dx  f  £gd£,  n  =  1(1)2 
0  0 


Reference  2  evaluates  these  thickness  integrals  for  a  typical  airfoil. 

From  the  oscillatory  coefficient  expressions  we  may  find  the  limiting 
values  and  the  steady  matrix  for  a  single  strip 

1(l+b/2d)  -b/d~j 
. -b/ 2d  b/dj 

[0  Bjfj]  f(l+b/2d)  -b/2dl 

0  DlfzJ  |_  "b/d  b/d  J 
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SECTION  II 


GENERAL  DESCRIPTION  OF  INPUT 


A.  Units 

Since  all  dimensional  input  is  geometrical  and  the  aerodynamic  matrix 
is  dimensionless,  only  a  consistent  set  of  length  units  is  necessary:  inches 
or  feet. 

B.  Classes  of  Numerical  Data  and  Limitations 

The  data  required  by  the  program  are  the  geometry,  reduced  velocities, 
and  Mach  numbers.  The  example  illustrates  their  use. 

1.  Example  Problem 

We  consider  the  four-strip  wing  shown  in  Fig.  2  at  the  freestream 
Mach  numbers  of  1.8  and  2.  5  for  the  reduced  velocities  of  4.  0  and  8.  0  and  the 
steady  case.  The  aerodynamic  matrices  will  be  found  with  and  without  the  tip 
correction,  and  the  thickness  integrals  will  be  generated  by  the  program  from 
the  expressions  of  App.  A  of  Ref.  2  for  an  assumed  airfoil  (constant  across 
the  span)  having  10  percent  thickness,  maximum  thickness  at  40  percent  chord, 
and  a  blunt  trailing  edge  that  has  a  trailing-edge-to-maximum-thickness 
ratio  of  0.  15. 

2.  Program  Restrictions  and  Options 

a.  The  number  of  strips  into  which  a  wing  may  be  subdivided 
must  be  <  25. 

b.  The  number  of  reduced  velocities  used  for  any  one  Mach 
number  must  be  <  50. 

c.  The  number  of  values  of  Mach  number  must  be  <  25. 

d.  If  it  is  desired  to  compute  the  steady  matrix  [C^g],  a 
zero  or  negative  value  of  V/brw  must  be  used  and  S  and  c  must  also 

be  supplied. 


1  1 


Strip 

No.  Ay  (ft) 

■b  (ft) 

d  (ft) 

1 

4.7 

12. 28125 

11.9 

2 

4.  2 

9.  50000 

9.0 

3 

3.  6 

7. 06250 

6.  6 

4 

3.  1 

4.  96875 

4.5 

b  =  ' 

6.  5  ft 

r 

s  = 

15.  6  ft 

c  = 

27.5  ft 

r 

Ct  = 

8.  0  ft 

tan  A  = 

0.  75  (cos  A  = 

0.80) 

tan  \  = 

-0.  50 

c  = 

2  1 . 0  ft 

S 

554.  0  sq  ft 
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Fig.  2.  Example  of  Four -Strip  Wing. 


Strip 

No. 

Av  (ft 

)  b  (ft) 

d  (ft) 

1 

4.7 

12. 28125 

11.9 

2 

4.  2 

9. 50000 

9.0 

3 

3.  6 

7. 06250 

6.6 

4 

3.  1 

4. 96875 

4.5 

b 

r 

6.  5  ft 

8  = 

15.  6  ft 

c  = 

r 

27.5  ft 

Ct  = 

8.  0  ft 

tan 

A  = 

0.  75  (cos  A  = 

0.80) 

tan 

X 

1 

o 

o 

c  = 

2 1 . 0  ft 

S  =  554.  0  sq  ft 
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SECTION  III 


DATA  DECK  SETUP 


A.  Loading  Order 

Input  decks  punched  from  keypunch  forms  are  loaded  behind  column 
binary  deck  HM10.  Any  number  of  complete  decks  may  be  stacked.  The 
data  for  each  deck  must  be  in  the  following  order: 

(1)  Heading  Card  1 

(2)  Heading  Card  2 

(3)  ISTRIP,  JMACHi  KVBRW,  ISTHK,  ISTIP,  NO  PUNJ 

(4)  b  ,  8,  S,  c,  c.,  tan  A  ,  tank,  cos  A 

r  t 

(5)  b^  series 

(6)  di  series 

(7)  Ay^  series 

(8)  M.  series 

J 

(9)  (V/bf  co)^  series 

(10)  Thickness  integrals 

(a)  Provide  thickness  integrals 

(b)  Input  data  to  compute  thickness  integrals 

B.  Input  Data  Description 

(1)  Heading  Card  1  may  contain  any  information  desired  (Columns  2 

through  72).  Column  1  should  be  blank. 

(2)  Heading  Card  2  provides  for  more  information  (Columns  2  through 
72):  vehicle,  date,  name.  Column  1  is  blank  and  all  or  any  of  the 
fields  may  be  blank  but  the  card  must  be  included. 
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(3)  Control  card  (FORMAT  1814) 


(a) 

1STRIP 

=  number 

of  strips 

(b) 

JMACH 

=  number 

of  Mach  numbers 

(c) 

KVBRW 

=  number 

of  reduced  velocities 

(d) 

ISTHK 

=  0  if  thickness  integrals  are  tabulated 

=  1  if  x  /c,  t,  and  r  are  constant 

m 

=  2  if  x  /c  and  r  are  constant  and  t 
m 

varies  with  each  strip 

=  3  if  x  /c,  t,  and  r  dll  vary  with  each  strip 
m 

(e)  ISTIP  =  0  if  tip  correction  factor  is  1.  0 

^  0  if  tip  correction  factor  i.3  to  be  computed 
by  program 

(f)  NO  PUNJ  =  0  or  blank  if  the  computed  matrices  are 

to  be  punched  in  cards 

i  0  if  no  punched  cards  are  desired 

(4)  Single  parameters  (FORMAT  6E12.  8):  When  computing  both 
steady- state  and  oscillatory  AICs  from  the  same  data  deck,  we 
use  the  eight  parameters:  br  (semichord),  s  (semispan),  S 
(surface  area),  c  (mean  aerodynamic  chord),  cj  (tip  chord), 
tan  A  (tangent  of  leading  edge  sweep  angle),  tan  X.  (tangent  of 
trailing  edge  sweep  angle),  and  cos  A  (cosine  of  leading  edge 
sweep  angle).  For  computing  oscillatory  cases,  S  and  c  may  be 
blank  (or  included);  and  for  steady  cases,  br  and  s  may  be  blank 
(or  included).  These  parameters  are  tabulated  in  the  same 
order  as  defined,  and  always  require  two  lines  (cards). 

(5) ,  (6), and  (7)  Series  bj,  dif  and  Ayt  (FORMAT  6E12.  8):  The 

number  of  lines  (cards)  required  by  the  b^,  d^,  and  Ay^  series 
is  controlled  by  the  number  of  strips.  Using  the  maximum  25 
strips  would  require  five  cards  for  each  series.  Tabulate  b., 
(local  semichords),  d^  (distance  between  control  points),  and 
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Ay^  (strip  widths)  in  the  order  defined  with  each  new  series 
starting  on  a  new  card. 

(8)  and  (9)  M^  (freestream  Mach  number)  and  V/b^ui  (reduced 

velocity)  series  (FORMAT  6E12.  8):  The  M.  series  may  be 

J 

input  in  any  order  desired  with  j  <  25  per  deck.  M.  cos  A 

J 

must  be  >  1.  25.  The  (V/b^cj)^  series  is  tabulated  only  once 
per  deck  and  will  be  repeated  by  the  program  for  each  Mach 
number  indicated  by  the  M^  series.  For  the  steady  case, 
tabulate  either  a  zero  or  negative  V/b^ox 

(10)  Thickness  integral  input  (FORMAT  6E12.  8) 


(a) 


(b) 


ISTHK  =  0;  the  thickness  integrals  must  be  tabulated 

in  the  order:  1^,  I^,  I I^,  Ig,  Jj,  J^,  Jg, 
J ,  K^,  K^,  Lj,  and  (three  cards). 


ISTHK  ^  0;  the  thickness  integrals  are  computed  by 

supplying  xm/c  (percent  chord  location  of 
maximum  thickness,  t  (airfoil  maximum 
thickness  ratio),  and  x  (ratio  of  Tte/T)* 

When  using  this  option  omit  (10.  a); 

1)  ISTHK  =  1;  the  three  ratios  are  each  constant  for 

the  surface  and  tabulated  in  the  order 
defined.  (First  field  in  three  separate 
cards. ) 

2)  ISTHK  =  2;  the  ratios  x  /c  and  r  are  each  constant 

m 

for  the  surface  but  t  varies  with  each  strip. 

Input  in  the  order  x  /c,  (Field  1,  first 

m 

thickness  integral  card),  t  series  (beginning 
on  a  new  card),  and  r  (Field  1  of  the  last  card). 
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3)  ISTHK  =  3;  the  ratios  x  /c,  t,  and  r  all  vary  with 

each  strip.  Tabulate  in  order  with  each 
new  series  starting  on  a  new  card. 

C .  Example  Keypunch  Forms 

Example  keypunch  forms  are  given  on  the  following  pages.  Columns 
73  through  80  are  reserved  for  data  deck  identification.  This  space  may  be 
used  in  any  fashion;  however,  it  is  suggested  that  the  last  three  columns  be 
used  for  sequencing.  Only  the  cards  with  sequencing  in  Columns  73  through 
80  are  to  be  used  in  the  sample  data  deck;  the  lines  (cards)  with  Columns  73 
through  80  blank  are  for  clarification  of  input. 
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80  COLUMN  FREE  KEYPUNCH  \At\ )  AEROSPACE  CORPORATION 
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18 


AEROSPACE  CORPORATION 
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80  COLUMN  FREE  KEYPUNCH  yfil  )  AEROSPACE  CORPORATION 


) 
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0  COLUMN  FREE  KEYPUNCH  J£\  AEROSPACE  CORPORATION 


SECTION  IV 


PROGRAM  OUTPUT 


Printed  Output 

1.  All  input  data. 

2.  Each  group  of  AICs  with  the  associated  Mach  number  and  V/b^w. 

3.  Sequencing  number  of  the  first  and  last  punched  cards  (output) 
for  each  group  (one  V/b^w)  of  influence  coefficients. 

4.  Example  problem  printed  output  is  given  on  the  following  pages. 


OSCILLATORY  CASE 
l/KIR)  »  0.80000000E  01 


» A 
00 
O' 
>0 
fNJ 

o 

oo 


o 

UJ 
O' 
« 
oo 
4 
rsi  f- 
m 
>-  rsi 
00  * 
• 

<M  O 

I 

It 

UJ 


UO  H-* 

rg 
«  O 


X  * 

Cl  Uf> 
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B. 


Punched  Output 


1.  A  deck  of  punched  cards  (output)  from  this  program  is  suitable 
as  an  input  deck  to  other  programs  requiring  the  use  of  AICs. 

2.  All  punched  output  is  sequenced  in  order  on  Columns  73  through 
80  starting  with  HM 100000.  The  data  appear  in  the  following 
order: 

a.  Card  1  contains  (V/b^o))^,  and  Mj  (FORMAT  6E12.  8). 

b.  Card  2  contains  m,  the  size  (number  of  control  points) 
of  the  AIC  matrix  and  n,  the  number  of  strips 
(partitions)  FORMAT  (1814). 

c.  The  AIC  matrix  punched  in  column  binary  form  and  its 
TRA  card  make  up  the  remainder  of  the  punched  output 
for  (V/b^co)^. 

3.  The  order  of  statement  (2)  is  repeated  for  all  Mach  numbers  and 
reduced  velocities  per  input  deck. 

4.  Each  matrix  is  punched  in  compact  form  by  columns.  Column  1 
begins  in  origin  1  and  Column  2  in  location  (1  +  matrix  size). 

5.  The  oscillatory  matrix  is  punched  in  the  order:  Column  1  (real), 
Column  1  (imaginary);  Column  2  (real),  Column  2  (imaginary); 
Column  m  (real),  Column  m  (imaginary).  In  the  steady  case,  all 
columns  are  real  and  are  punched  in  order. 
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SECTION  V 


PROCESSING  INFORMATION 


A. 


B. 


Operation 

STANDARD  FORTRAN  MONITOR  system 
Estimated  Machine  Time 
T  =  time  In  minutes 

ISTRIP  =  number  of  strips 

JMACH  =  number  of  Mach  numbers 

KVBRW  =  number  of  reduced  velocities 

n  =  number  of  sets  (decks)  of  input  data 


T  =  .5  +  ..  03  [(ISTRIP  ■  JMACH-  KVBRW  ^  +  (ISTRIP  •  JMACH  •  KVBRW)2 
■+  •  +  •  +  (ISTRIP-  JMACH-  KVBRW )J 


C.  Machine  Components  Used 
Core  storage:  6772 

Standard  FORTRAN  input  tape  (NTAPE  2) 
Standard  FORTRAN  output  print  tape  (NTAPE  3) 
Standard  FORTRAN  punch  tape  (NTAPE  7) 
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SECTION  VI 
PROGRAM  NOTES 


A.  Subroutines  Used 

RDLN:  reads  and  prints  title  cards 
AEROP3:  punch  AIC  matrix 
BINPU:  column  binary  punch 

All  other  subroutines  used  are  on  library  tapes. 

B.  Generalized  Tapes 

Input,  print,  and  punch  tapes  in  this  coding  are  defined  as  Units  2,  3, 
and  12  respectively;  however,  these  may  be  altered  by  placing  the  desired 
units  on  symbolic  cards  HM100054,  HM100055,  and  HM100056. 


SECTION  VII 


FLOW  DIAGRAM 


A  partial  list  of  the  principal  FO 
may  be  related  to  the  physical  notation 


FORTRAN  Notation 

BI 

DI 

DEFY 

XMACHS 

VBRW 

XMC 

TAU 

AR 

DISTO 

DISTI 

CHI 

XLHO 

XLAO 

XMAO 

XMHO 

All,  AI2 ,  .  .  .  ,  AIN 
A  J 1 ,  AJ2,  .  .  .  ,  AJN 


AK1 ,  AK2,  .  .  .  ,  AKN 


VIII 
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SYMBOLIC  LISTING  (continued) 


FORTRAN  Notation 


Physical  Notation 


AL1 ,  AL2,  .  .  .  ,  ALN 
B 

BB1,  BB2,  .  .  .  ,  BBN 
Cl,  C2,  ...»  CN 
Dl,  D2,  .  .  .  ,  DN 
FF1 ,  FF2,  ....  FFN 
ISTRIP 
JMACH. 

KVBRW 

ISTHK 

IS  TIP 

NO  PUNJ 

BR 

SMS 

CAPS 

CBAR 

CT 

TAN  LA 
TANLAM 
COS  LA 


LP  L2’  *  •  •  ’  Ln 

P 

BL  B2 ’  •  '  V  Bn 
CP  C2’  •  '  '  '  Cn 

Dp  D2"  ‘  ’  *  Dn 

f  f  f 

P  2* -  '  ’  n 

Number  of  strips 
Number  of  Mach  numbers 
Number  of  reduced  frequencies 
Thickness  integral  control 
Tip  correction  control 
Punch  control 
b 

r 

s 

S 

c 

Ct 

tan  A 
tan  X. 
cos  A 


The  complete  symbolic  listing  is  given  on  the  following  pages. 
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93  WRITE  OUTPUT  TAPE  NTAPE3*  16  HM100493 
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